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Abstract 

Within the framework of simple perturbation theory, recurrence time of quan- 
tum fidelity is related to the period of the classical motion. This indicates 
the possibility of recurrence in near integrable systems. We have studied such 
possibility in detail with the kicked rotor as an example. In accordance with 
the correspondence principle, recurrence is observed when the underlying clas- 
sical dynamics is well approximated by the harmonic oscillator. Quantum re- 
vivals of fidelity is noted in the interior of resonances, while classical-quantum 
correspondence of fidelity is seen to be very short for states initially in the 
rotational KAM region. 

I. INTRODUCTION 

Qualitative dynamical behaviour of a classical system is characterized by the sensitivity 
of its orbits with respect to initial conditions. The sensitivity also changes with parameters 
of the system. If an orbit is not sensitive to initial conditions, the underlying dynamics 
is regular. On the other hand, for chaotic dynamics the orbit is highly sensitive to initial 
condition. Lyapunov exponents, which measure the rate of exponential divergence of two 
initially close orbits, is zero for the regular case and positive for the chaotic case. An 
analogous classification of dynamics is not possible in the quantum domain, since quantum 
theory does not accommodate the very concept of orbits. 

If we associate two localized wave packets which are initially close-by in state space, over- 
lap of these states remain invariant under unitary quantum evolution. The corresponding 
classical Liouville evolution of phase space densities is also a linear unitary one, but the abil- 
ity to develop structures on infinitesimal scales, without the restrictions of the uncertainty 
principle, coexists with underlying orbit chaos. In the quantum mechanical case exponential 
instability of for instance the wave packet center is seen for quantized chaotic systems over 
a short time scale, corresponding to the Ehrenfest time scale. 
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One approach to quantify quantum sensitivity was initiated a while back by Peres [1] 
who proposed the overlap of states that are evolved from the same initial state \a) under 
two Hamiltonians differing by a small perturbation. The overlap intensity, also known as 
fidelity, is defined as 

f(t) = \{a\e iHlt/n e- m>t/n \a}\ 2 (1) 

where Hq represents the Hamiltonian of an unperturbed system, H\ = Ho+eV represents the 
perturbed system and e is the perturbation parameter. Fidelity has evoked enormous interest 
in recent years with different interpretations. It is a measure of quantum irreversibility in 
the context of decoherence when the system under investigation interacts with a classically 
chaotic environment [2]. It is realized as polarization echo in nuclear magnetic resonance, 
measuring the irreversibility in many body quantum dynamics [3]. It also characterizes the 
loss of phase coherence in quantum computation [4]. 

Initial investigations by Peres [1] show that f(t) is appreciable (on a time average) and 
its fluctuations are more for the regular case than for the chaotic case. Some of the recent 
studies have focused on the decay rate of f(t) for chaotic systems [5]. More investigations 
on this topic show that while the fidelity decays on time scale ~ 1/e 2 for chaotic systems, it 
decays on shorter time scale ~ 1/e for integrable systems [6] . In [7] the decay for integrable 
system is shown to have power-law behaviour. 

In this paper, using perturbative approximations, we relate recurrence time of the quan- 
tum fidelity to the classical period. This indicates that recurrence of quantum fidelity would 
be possible in nearly integrable systems. With the kicked rotor as a dynamical model, we 
examine this recurrence phenomena in detail. Our numerical experiments also reveal that 
the evolution of fidelity exhibits a wide range of behaviour, depends on the choice of initial 
state. 



II. FORMULATION 

Let us consider an unperturbed system H = H(u), with uj being the system parameter, 
satisfying the eigenvalue equation 

H \u n ) = E n \u n ) . (2) 

Let Hi = H{uj') be the perturbed system where u' = uj + e and e is the perturbation 
parameter. Under a small perturbation, for the simplest zeroth order approximation to the 
fidelity it is sufficient to assume that the energy eigenvalues differ up to the first order of e 
and the corresponding eigenstates remain unaffected. It is easy to see that the eigenstate 
changes contribute to the fidelity at an order e. 

That is the approximation we can make for the calculation of fidelity amounts to the 
approximation: 

dE n 

« E n \u n ) ; E n ^E n + e—^. (3) 

aw 

Here dE n /dui measures the rate at which nth level changes with the parameter. It is also 
known in the literature as level velocity wherein the parameter is thought of as pseudo-time. 
In this approximation the fidelity can be written as 
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where c n = (u n \a). Note that /(0) = 1 and we are interested in the time evolution of the 
fidelity. Assuming that the weighting probabilities \c n \ 2 are strongly centered around a mean 
n, one can expand dE n j duj using Taylor series in n around n up to the first-order correction 

as 



dE n 



cku duj 
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(5) 



This gives the recurrence time 



d?E„ 



dnduj 



(6) 



such that the fidelity becomes unity when t is an integer multiples of t r . Recollecting that 
the classical action S is quantized as (2nh)n = S, 



d 2 E n 



dS dco 



(7) 



In terms of T = dS/dE, the period of the classical orbit underlying the coherent state \a), 
the recurrence time is thus given by a combination of simple perturbative and semiclassical 
approximations as 



e 



dT 



duo 



(8) 



Thus the recurrence time of quantum fidelity is directly related to period of the underlying 
classical motion as well as to its variability to the relevant parameters. 



III. HARMONIC OSCILLATOR 

According to the above relation, the quantum fidelity for the harmonic oscillator 

H(u) = t. + i w y (9) 

with a classical period T = 2tt/lj recurs at 

*r = Y ■ ( 10 ) 

This recurrence time can also be obtained from the underlying classical dynamics by evolving 
an initial phase space density p (q, p) which satisfies the normalization condition 

J p(q,P)dqdp= 1. (11) 



3 



Analogous to the quantum fidelity, we define normalized classical fidelity as 



F(t) = jy p 2 (q, p) dq dp J J p (q, p) p (q 2t , P2t) dq dp . 



(12) 



where p (q2t,P2t) is the final phase space density which is obtained from forward time evolu- 
tion of the initial density for a time t with the frequency being u, followed by backward time 
evolution for the same time t with the frequency set as a/. The classical fidelity measures 
the overlap of the initial density with the final density so obtained. 

If the forward time evolution of the oscillator is given by (q,p) — > (g, p) i.e., 
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(13) 



the initial density evolves forward in time t as p (q,p) — > p (qt,Pt) where q t and p t are given 
by the relation: 
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(14) 



Then we consider the backward evolution for the same time t with the frequency now set to 
P (Qt,Pt) — >• P (q2t,P2t) where q 2t and p 2t are given by the relation: 



q2t 




_P2t _ 





cos(o/t) (l/o; / )sm(a; / 
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We shall now write 
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(16) 



where the elements of square matrix are obtained by combining the Eqns. (14) and (15). 

If q2t = q and p 2 t = P for some time t = t^, then F(tn) = 1 and the classical fidelity 
recurs fully. This condition is given by 



a-l b 
c d-l 



0. 



(17) 



Since (ad — be) = 1, the above condition becomes a + d — 2. The recurrence time t R then 
satisfies the relation 



2 cos(o;t_R) cos(u't R ) + 



u 2 + uj' 2 



UJLd' 



sin(o;tft) sm^u'tn) = 2 



(18) 



Substituting u' = u + e and approximating u' 2 ~ c<j 2 + 2ec<j (up to the first order of e), this 
condition reduces to 



/ x , 27r 

cos(et R ) = 1 or t R = — . 



(19) 
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Thus recurrence time t R of the classical fidelity, at this level of approximations, is identical 
to that of the quantum fidelity given by the Eqn.(lO). This shows that recurrence of quantum 
fidelity for the harmonic oscillator is a manifestation of the underlying classical dynamics. 
If we replace u by iuj and ui' by iu>', the dynamics corresponds to the unstable motion of an 
inverted harmonic oscillator. In this case the recurrence condition is 



2 cosh(ut R ) cosh.(u't R ) — 



LO 



UJUJ 1 



sinh(u;tR) sinh(a/ifl) = 2 . 



(20) 



Up to the first order of e, this condition reduces to cosh(et^) = 1 which has only the 
trivial solution tR — 0. Hence there is no recurrence in the classical fidelity of the unstable 
oscillator. 



IV. KICKED ROTOR 

For further investigations on the recurrence phenomena in near integrable systems, we 
consider the kicked rotor whose Hamiltonian is 

H = £ + V(q)526U-t). (21) 
A i 

Here V(q) = k cos(2nq) / 4n 2 and k is the kick strength - the only parameter of the system. 
The corresponding kick-to-kick dynamics is given by the standard map: 

Pt+i = Pt + (fc/27r) sin(27r % ) \ 

> (mod 1) . (22) 
qt+i = qt + Pt+i J 

Note that the modulo operation restricts q and p to between —1/2 and 1/2 with the opposite 
edges being identified. The kick-to-kick quantum dynamics is given by the quantum map: 

|V t+ i) = U \A) I U = exp[-tp 2 /2h] exp[-iV(q)/h] (23) 

where Uo = U (k) is the unit-time quantum propagator. 

The 2-torus phase space can be quantized upon introducing periodic boundary conditions 
in both the canonical variables [8]. This imposes finite number, N, of quantum states such 
that iV = l/2n% (N — > oo is the classical limit); we take iV = 500 for the following calcula- 
tions. Since the Hamiltonian is time periodic, according to Floquet theory [9] fundamental 
solutions statisfy the eigenvalue equation 

E/o|0n> = e*»|0„) (24) 

where n are real (between and 2n). If n = E n /h, then E n are called quasienergies which 
are analogous to the energy eigenvalues of the time-independent Hamiltonian system [10]. 
Correspondingly the states \<p n ) are called quasienergy states. It is worth noting that as a 
consequence of hermiticity of the Hamiltonian, the quasienergy states \<p n ) are orthogonal 
and they form a complete set in the N dimensional space. The quasienergy spectrum can be 
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numerically obtained by diagonalizing the matrix representation of U. Choosing the discrete 
position eigenstates |Z) as the basis, [/-matrix takes the form [11] 



(l\U\l') = -^exp 



(25) 



where /, V = -N/2, -N/2 + 1, . . .N/2-1. 

Let us consider U as an unperturbed system and U\ = U(k + 5k) as the perturbed 
system. Then the quantum fidelity is given by 



/(f) = K«|tfr*W 2 



(26) 



In accordance with our earlier approximations the eigenstates of the perturbed system are 
taken to be such that 



£/i|<M~e^|0 n > ; 



+ 5k- 



dk 



(27) 



where d(j) n /dk = h 1 ((j) n \dV/dk\(j) n ), quantum fidelity for the kicked rotor can be approxi- 
mated as 



/(*) 



J2 |(0«|«)| 2 exp 



it5k 



h \ dk 



(28) 



This gives the recurrence time as 



5k 



dT 



dk 



(29) 



with T being the classical period of the rotor. 



V. NUMERICAL RESULTS 

The dynamical features of the standard map depends on the kick strength. In the absence 
of kick (k = 0), it is just the twist map with regular dynamics. For small kick strengths 
(k < 1), the phase space is predominantly regular with primary nonlinear resonance zones 
along p = and there are large number of smooth Kolmogorov-Arnold-Moser (KAM) tori 
in large \p\ regions on which the motion is quasiperiodic (see Fig. 1). These two distinct 
phase space regions are demarcated by a separatrix. In this case, the phase space is similar 
to that of the simple pendulum. As k increases, the KAM tori are slowly destroyed and 
at k ~ 1 all the KAM tori are destroyed [12] leading to the onset of global chaos. When 
k ^> 1, the dynamics is predominantly chaotic. In this paper, we take k = 0.3 and 5k = 0.01 
throughout, as our interest is restricted to the near integrable regime. 

To evaluate the quantum fidelity for different regions of the phase space, we take the 
initial state as a coherent state peaked at (qo,Po), i-e., \a) = \qo,Po). It is a minimum 
uncertainty (very nearly Gaussian) wave packet with equal spread a = \Jh/2 in both the 
canonical variables. Such an initial state on the torus can be constructed by the method 
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devised by Saraceno [13]. The initial states we have chosen for the following computations 
are marked in Fig. 1. 

For comparison, we also compute the classical fidelity by taking an initial phase space 
density to be a Gaussian 



2a 2 



\ {(q - ?o) 2 + (p - Po) 2 } 



(30) 



This density is a classical equivalent of the coherence state \qo,Po) containing an ensemble 
of 10 4 initial conditions. For the classical evolution, each point of the initial density is 
forward iterated using the standard map (with parameter k) for t steps, followed by backward 
iteration for t steps by the corresponding time reversed map, with k being replaced by k + Sk. 
This gives the final density p (q2t,P2t)- With this, we use a discretized version of Eqn. (12) 
to calculate the classical fidelity. 

It should be emphasized that the recurrence of fidelity is possible only when the pertur- 
bative approximation is valid. In order to validate this approximation for different regions of 
the phase space, we compute phase space representation of inverse participation ratio (IPR) 
as 

^,p)=£|(0nk,p)| 4 (31) 
n 

such that z~ x measures the effective number of states \(f> n ) required to construct a given 
region of phase space. As we see from Fig. 2, while very few states are required to construct 
the region near the stable fixed point, many states contribute to construct the separatrix 
region and rotational KAM region. In other words, low energetic part of the spectrum 
that corresponds to the neighborhood of stable fixed point has very low density of states; 
high energetic part of the spectrum that corresponds to the separatrix and KAM regions 
are highly dense. Generally, such simple perturbative expansions work badly for the dense 
spectrum. Hence the approximation is expected to be valid near the stable fixed point, and 
so is the recurrence of quantum fidelity. For other cases, the approximation would not be 
sufficient and the recurrence is not expected. As we demonstrate below, this is indeed the 
case. Nevertheless, for the dense spectrum, "quantum revivals" of fidelity may occur which 
is beyond the scope of present study, but it is under active study. 

To begin with, we focus on the neighborhood of the stable fixed point. Note that Jacobian 
of the standard map evaluated at (0.5, 0) has eigenvalues e ±tX and A is obtained from the 
relation 

e iX + e~ iX = 2 cos X = 2 — k . (32) 
Since A is the frequency of the oscillatory motion in this region, the corresponding period is 

-i 



T(k) = ^ = 2n 
A 



cos 



(33) 



and we obtain the recurrence time as 



2tt 



t r = ^k(4-k). (34) 
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For the parameters we have considered here (k = 0.3, 5k = 0.01), the quantum fidelity recurs 
at t r = 662. If we approximate the underlying oscillatory motions by harmonic oscillators 
with frequencies 



cos — - — 



uj' = oj + e 



cos 



(35) 



then t r = 2n/e = 667. 

In Fig. 3 we show the evolution of fidelity for the initial state |0.45,0) that is localized 
in the neighborhood of the stable fixed point. For initial time steps, quantum fidelity has 
a smooth fall-off from unity to zero. It remains nearly zero for some time and then it rises 
and recurs at t = 667. This is close to the value 662 that is obtained from perturbation 
theory. Notice that the recurrence time matches exactly with that of the harmonic oscillator 
approximation. The fidelity also recurs periodically and the agreement with perturbation 
theory is found to be good. In addition, we find that the classical evolution exactly retraces 
quantum evolution. Thus, fidelity for oscillatory motion in the neighborhood of a stable 
fixed point exhibits very good quantum-classical correspondence. 

Fig. 4 corresponds to the state which is initially placed in the interior of the primary 
resonance zone. In this case, the quantum fidelity almost recurs i.e, 95% recurrence is 
observed at t = 670. Note that this time is close to the recurrence time predicted from 
the harmonic oscillator approximation. Here also the perturbative approximation agrees 
fairly well with the actual evolution. However, unlike the previous case, classical evolution 
initially follows the quantum evolution and then deviates at later time steps. This shows 
that the observed near recurrence of fidelity is a "pure" quantum phenomena. This shows 
interestingly that, within the primary resonance zone, the quantum dynamics seems to be 
less sensitive to nonlinearity. 

As we see from Fig. 5, quantum evolution is highly complex when the underlying classical 
motion is on the separatrix. The perturbative approximation fails and the fidelity does not 
recur. In this case, asymptotically the classical fidelity becomes more stable in comparison to 
the quantum counterpart. If we approximate unstable motion on separatrix by an inverted 
harmonic oscillator, classical fidelity also does not exhibit any recurrence. 

Finally, we move on to the rotational KAM region of the phase space. A typical evolution 
of an initial state localized in this region is shown in Fig. 6. As expected, perturbative 
approximation fails and there are no signatures of fidelity recurrence. Although classical 
motion on the KAM tori is regular, we notice that quantum-classical correspondence is seen 
only for a short time. This may be due to enhanced quantum interference in this region of 
energetic rotational motion. More study on these issues is currently underway. 



VI. CONCLUSION 

We have shown that, within the framework of simple perturbative and semiclassical 
approximations, recurrence time of quantum fidelity is related to period of the classical 
motion. In the case of the harmonic oscillator, fidelity recurs at t r = 2-rr/e where e is the 
perturbation parameter. Invoking the kicked rotor as an example, we have studied in detail 
the recurrence phenomena for nearly integrable regimes. It is shown that quantum fidelity 
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recurs only when the underlying classical motion is well approximated by the harmonic 
oscillator. This recurrence is also in accordance with the correspondence principle. Though 
the recurrence is not possible in other cases, quantum revival of fidelity can occur and this 
needs further investigation. 
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FIGURES 
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FIG. 1. Phase space portrait of the standard map for k = 0.3. The points 
(lOiPo) = (0.45, 0), (0.3, 0), (0, 0) and (0,0.3), shown as solid circles, correspond to the intitial 
coherent states used for the quantum evolution. 
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FIG. 2. The phase space representation of the Inverse Participation Ratio z(q,p) plotted for 
k = 0.3. 
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FIG. 3. The solid line is the quantum fidelity for the initial state centered at (0.45, 0), near the 
stable fixed point, and the dashed line is the perturbative approximation (28). The corresponding 
classical fidelity, evaluated from Eqn. (12), is shown as dots. 
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FIG. 4. The solid line is the quantum fidelity for the initial state centered at (0.3,0), in the 
interior of the primary resonance zone, and the dashed line is the perturbative approximation 
(28). For comparison, the corresponding classical fidelity is shown as dots. Inset shows the initial 
evolution. 
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FIG. 5. The solid line is the quantum fidelity for the initial state centered at (0, 0), the unstable 
fixed point, and dashed line is the perturbative approximation (28). Dots are the corresponding 
classical fidelity and the initial evolution is shown in the inset. 
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FIG. 6. The solid line is the quantum fidelity for the initial state centered at (0,0.3), on a 
rotational KAM curve, and the dashed line is the perturbative approximation (28). Dots are the 
corresponding classical fidelity and the initial evolution is shown in the inset. 
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